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ABSTRACT: 

A study of the effects of using mean absolute deviation (MAD) 
to estimate variability in setting reorder levels f"*!* the inventory 
of a stock item. The method presently employed by NavSup in setting 
such reorder levels involves exponentially smoothed estimates of 
the mean and variance of the demand process. Any error Involved 
in setting reorder levels results in a change in the underlying 
risk which in turn can be translated into costs. Such errors for 
the method of estimation presently employed are compared with 
standard maximum likelihood procedures. By simulating several 
normal systems, the smoothing technique is found to be inferior to 
classical methods with no reduction in computational difficulties. 
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1. INTRODUCTION AND SUMMARY 

In a previous report the use of MAD (mean absolute deviation) 

in the analysis of demand variability as proposed by R. 6. Brown [3] 
and presently used by NavSup [4] was discussed. In particular, the use 
of exponentially smoothed estimates of MAD was criticized on theoretical 
grounds and a case was made for further study Into the probability 
distributions of the various estioiators presently in use. Of special 
importance, the use of MAO in setting safety levels was mentioned as a 
problem area in need of further investigation. 

The present report concerns Itself with the extension of research 
Investigations in the areas just mentioned. Several attempts to derive 
the necessary probability distribution theory for MAD estimators 
resulted in disappointing futility for the most part. This was 
anticipated for, as previously mentioned in [1], absolute deviation as 
a measure of variability has been abandoned for years on theoretical 
grounds by Statisticians. 

Two minor results that were obtained concern the ratio of MAD, A, 
to standard deviation, a. The complete result of this ratio has been 
determined for the Poisson case. As alluded to in [l], the universal 
use of ^ 0.8 for this ratio is especially bad for low demand 

Poisson item. The same general criticism can be made in the negative 
binomial case. However, it was found that for the choice of parameters 
presently being used by NavSup, the approximation of this ratio by 
0.8, is a good one. 
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Falling to obtain the required distribution theory to study the 
behavior of HAD estimates, simulation was used to compare smoothed 
estimates of o to the usual maximum likelihood estimates mostly for 
the case of normal demand. Parameter choices that reflect high demand 
Items were selected in order to give MAlB the best possible advantage In 
the comparison with alternate methods of estimation. Such parameter 
choices were made after consultation with NavSup personnel In order to 
duplicate situations that actually exists moreover, large amounts of real 
data are available for such cases. In every single case studied, the 
sample variance of the smoothed estimate of MAD and consequent estimate 
of cr was roughly twice that of the maximum likelihood estimates of (X. 
The same was true In a slow mover case or two that was examined. This 
confirms In a general way the observation made in a similar study carried 
out by Asher and Wallace C2] In which they found HAD to be significantly 
less efficient (about 20^ compared to the classical minimum variance 
estimators under the 6auss~Markov assumptions. 

Not content merely to confirm this observation, however, the estimates 
found were used further to see what the actual effects on reorder levels 
were* This was approached in several ways, first, MB and maximum 
likelihood were compared by computing the percentage of tim© the known 
theoretical reorder level was overestimated, thereby resulting in too 
much on hand stock. Secondly, the results were examined to see what 
percentage of the observations fell within k units of the theoretical 
order level for k ^ 1,2,*** and various fixed values of risk. Finally, 
the two methods were compared by setting reorder levels and then computing 



the actual risk attained by those levels where.j, for each case agaln> the 
theoretical or true risk Is known. As seen In Section 4^ exponentially 
smoothed mean and MAD came off second best compared with maximum likeli- 
hood estimates of M- and c in every single case examined . 

Recognizing the limitacions of simulation, no sweeping claim is here 
made for proof that exponentially smoothed MAD is an Inefficient method 
of accounting for demand variability. At the same timei, one cannot ignore 
the fact that present methods were uniformly Inferior In the situations 
examined . 

In Section 5 recommendations for further study are made. Among these 
is the suggestion that real data be used from the histories available at 
FMSO for comparing the results in retrospect with what would have been the 
case had maximum likelihood procedures been used. The writer wishes to 

oo 

thank Lt. Ozden Orneck for his Invaluable assistance in constructing and 
running the computer programs as well as assistance In mathematical 
derivations. Acknowledgement should also be given to Cdr. Jack E. White of 
FMSO for his unfailing cooperation in defining the problem areas and supply- 
ing parameter values that are realistic in terms of NavSup use. Credit 
should also be given to Mr. James W. Prichard (SUP 04E) for his 
continued endorsement and interest in this research area. 

2. RATIO OF A TO g 

It will be helpful to review seme parametric definitions and 
establish a notation to be used here and in ensuing sections. Let X 
be a random variable with mean M> and variance (standard deviation 

o) o The mean absolute deviation, MAD for short, will be denoted A 
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and is defined as 



(2.1) L * ECjl'nil) 

As remarked in [iJs It w©wld be m®re rational to define MD as E(|x-in|) 
where m is the median of X bait It is (2.1) tha.t is ®sed by Brown [3] 
and NavSup [4]. Gonse^iuently the same definition is adopted in this 
report. 

Now when X is normally distribnteda it is well known that the 
ratio of A to 9 is given by /t roughly 0,8. It is somewhat 

surprising that the ratio is approximately the same for certain other 
families of probability distributions such as the Exponential!, Uniform 
and Triangular. Howevers, this result is not universally true and is 
particularly a poor approximation for the Poisson family j, a model often 
used for the so°called slow mover type of inventory item. This fact 
was demonstrated earlier in [l3 by examining the ratio for selected 
values of the Poisson parameters X, That ratio has now been determined 
for all values of X and it nay be instructive to see the behavior- of 
this ratio in a complete sense. 

By [z], we shall mean the greatest integer in Sj, that is, the 
largest integer n such that n ^ a (and hence z < n + 1} . We 
denote the Poisson (X) mass function by p so that, 
p(x> » 

Then we have, recalling that p, ® “ X for this case, 

® [X^ *. 

A ® B(|x*=p 1]) * V |x°X|p(x]) ® ) CX=x)p(x) + y (x“X)p(x) 
x^ x-0 k=[x3+1 

z 

But, letting F(z) “ L p(x) define the Poisson distribution function. 



CX] CjO 

X (^“x)p(x) “ ^ F([x]) - ) xp(x> « X F([X] 



09 

X (X“^)p(*) “ y *P(^ - X + X F([X]) 
x*CX]+l x“[X]+l 



Adding these results. 






(2.2) A • 2 X F([X]) - 2 X + 2 A xp(x) 

x*CX]+l 



+ ) xp(x) - X 

X“[Xl+l 



But, for X 2 1, 

^ V X* 

xp(x) = xe A 

Hence, 




(x-l)i 



Xp(x-l) . 




xp(x) 

+1 



m 

^ T p(**^) 

x“C\3+x 



X y p(y) = X(1-F(CX]-1». 
y=tx] 



Substituting this result in (2.2), 

A « 2 X F([X]) - 2X + 2X-2X F([X]-1) 
= 2 X(fCX]-F(CX]-1)) 

“ 2 X p(CX]) 



Thus, 

(2.3) 



-X vW 

A “ 2 Xe X for any X > 0. 

rw 



Since o = /X we have, 



(2.4) 



^ = 2/X 



if 

Cx]j 



The graph of ~ as a function of X appears in Figure 1. 
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Figure 1. Graph of ^ 

The aitalysis that validates the graph is the following. 

If n<X<n+l then CX] “ n and, ^ x”^ e . 

w n • 

As a function, say f , of X in this open Interval, f is 

n n 

differentiable with, 

f'(X) X“"^ (n-l^-X) 

n 

Hence, f ^(X) >0 if X<n + % and f^(X) <0 if X>n + % 
n n 

so that f^ possesses a relative maxiinum at X n + % with a 

maximum value of fj^(n-t%)^ for n = 0,1,2,***. 

With the formulas established for f and f it is easy 

n n+1 

A 

to verify that ~ is continuous at each of the integers although 

not differentiable there. In any case, it is shown in the appendix 

, A 

that the sequence of values of ~ at n. 



n4% 

n 



-n 

e 
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is monotone increasing. Using Stirlings approximation [7] to nl , 



l2n-r 

we see that, 



/i — h — < f (n) < /i 

IT Hr TT 

I2n-1 

and lim £ (n> /2 . Thus the values of ~ converge in the limit 

n-»» “ n 

to the factor . However, as previously remarked in [l] and as 
TT 

seen visvalfy in Figure 1, approximating k. hy /i is particularly 

^ TT 

bad when X < 1 which is precisely when the Poisson assumption is 
of primary interest. 

For intermediate movers, the negative binomial distribution is 
employed as a model for demand. Specifically, from [4], such an 
assumption is made whenever mean demand is between 2 and 20 units. 
The negative bincmiial distribution is a family of probability 

■* i 

mass functions having two parameters r and p where r is a 
positive integer and 0 < p < 1 . For our purposes the general 
formula for the mass function is defined by, 

P(x) = pV . X “ 0,1,2,“* (q=l-p) 

The mean ^ and the variance of a random variable having this 
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mass function are given by. 



£a 



13 . 



Proceeding as in the Poisson case, 

Cii] 



A ■ E(|x-t*|) “ X Ix-^iIp(x) - V(p.-x)p(x) + y (x“p.)p(x) 
x^ x^ x=[flQ+l 



But, 



Cli] 



Cp.] 



I' + I (ii-*)('l;-Sp'< 



x-O 



X =1 



and. 



[p: 

I (X-P) ('-^-SpV ■ - I x('^-1)pV + P I 

.fij+1 x»l * x=0 



x*Cpt 

Adding these results. 



[pt] [^] 

(2.5) A = aCp.p’^ + M- ') p(x) - X xp(x)3 



xVl 



x«l 



_ , it is no problem to cocnpute the ratio -* for various 
P a 



expressing A in terms of finite sums. 

Since o = 

choices of r and p on the computer. 

It is easy to choose values of r and p for which the ratio 

is very different from 0.8. However, the interest would only be 

academic for values of r and p such that ^ ^ [2,20]. For that 

range of values, a total of 1,000 cases were considered as follows. 

Let \i € [2,20] be selected and for each such |Ji, r was allowed 

to vary over the set [2,4,6, *** ,20} . For each such choice, p is 

r 

then determined from the relation p and q is specified by 

r+n 



« 8 " 



was 



q “ 1 - p. Then A was computed using (2.5) and finally ^ 
evaluated* In each such **group** of ten parameter pairs, it was 
observed that ^ was invariably monotone increasing although there 
was no monoticity from group to group* More significantly however, 
it was observed that ^ varied always between a minimum value of 
0.750 and a maximum value of 0.805. For example, if p = 0.893 
and r “ 20 whence p. “ 2.4, we find that ^ = 0.805 although for 

P = 0.455 and r = 2, resulting also in n = 2.4, ^ » 0.770. On 

the other hand, for p = 0.174 and r * 4, yielding n = 19, 

^ •* 0.781 so that low values of ^ are not always associated with 

low values of pi. 

In summary, then, we may say that for the negative binomial 

case and the range of parameter values presently used by NavSup, 

the approximation /l for ^ appears to be a safe one. 

tT ° 

3. ESTIMATING VARIABILITY WITH MAD 

Having established the value of ^ for cases of particular 
interest in the inventory models, it follows that a = k A for 
some constant k. Consequently, if A is estimated directly from 

tv 

the data to obtain an estimate A , say, then a can be estimated 
by applying the formula a « k A . In this way, one obtains a 
measure of the variability of demand through estimates of A . 

How good such a procedure might be clearly depends upon the 
technique used to estimate A in the first place and what the 

tv 

effect is on scaling that estimate by the factor k to obtain a. 
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There are of course many different ways to estimate a . 

Following Brown [33 i> NavSup uses a formula based upon exponentially 
smoothing certain error forecasts defined as follows. Let Xq,Xj^, * * * ,X^ 
denote Independent demands over discrete units of time up to and 
Including t. Suppose that is a forecast of the demand at 

time t, such forecast being made at time t**l. Then the difference 
e^ = Xj. • ^«i 9 represents a forecast error that Is dependent upon 
the method of forecasting used. Again following Brown, suppose that 
exponential smoothing Is used to forecast so that, using smoothing 
constant or, 

t°2 



In this report we will only be considering the model where random 

2 

demand X is normal with constant mean ^ and variance cr . With 
this In mind we have, as previously reported in [l3, 

B (e^) = 0 for all t 

and the error terms have time°dependent variances given by. 









,crt-2p' 



2t»l 



However, if we let t ^ ® , we obtain a limiting variance of 



(3-1) 




e 




or 



a 



_ 

2 



a 

e 



Now for each t, e^, being a linear combination of Independent 
normal random variables. Is again normal and Its mean absolute 
deviation, Ag, Is related to by = ''■2 Ag as previously 



- 10 - 



established. Consequently, from (3.1) we may write, asymptetically. 



(3.2) 




This formula relates demand variability to the mean absolute 
deviation of error forecasts. 

The formula used in [3] and [4] to estimate 4^ and ultimately 
o is based on an exponentially smoothed function of the error 
forecasts. If 



then a is estimated for demands by 



As remarked in [l], the worth of this estimation procedure is 
difficult to assess without some knowledge of the probability 
distribtulon of . Serious efforts notwithstanding, that 
probability distribution has not been derived to date. The real 
basic difficulty is that random variables 
compose A^ , while normally distributed, are not mutually 



independent. Even the matter of finding the second moment of A^ 
has been Intractible thus far. 

In an effort to discover some Indication of the variability 
involved, a simulation study was undertaken and resulted in mounting 
evidence that there is reason to believe that the variance of A^ 
is bounded below by a positive quantity so that even in the limit 



t“l 




11 => 



the actual estimates o based on this quantity will fluctuate 
about the true value a. 

In turn, this immediately suggests comparing this method of 
estimation with the classical maximum likelihood estimate cr 
given by the formula. 



(3.4) 



a 



= i=l *■ 



2 



t 



Now, properties of a are well known (see [s] for example) , being 
derived from the so-called Chi-distribution. In fact. 



E(a) 







where 



rj- rx|) j 

2 



and. 



V(a) = ^ ^2 ^ where = 1 - 

t t t t 2(t-l) 

From these formulas we see that 

lim E(ff) = a 
t=»® 

(3.5) 

lim V(o’) = 0 

t-*« 

Hence, if it is true that the variance of o is bounded away from 

A 

zero, then, even though unbiased in the limit, the fact that o has 
a vanishing variance would make it preferable as an estimating tool. 
Our simulation results certainly seem to concur that this is the 
case. 

For simulation purposes, several members of the normal family 
were selected so as to be representative of a wide class of fast 
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moving Inventory Items. Using random generation samples of size 
200 were drawn from these populations and a running account was 

<K 

kept of the estimates o and O of the triie and known standard 
deviation < 7 . This was done by computing the values of (3.3) and 
(3.4) at each time period t. 

Incidentally, it should be observed In this regard, that the 
running computation of 9 involves no more time nor storage than 
that of 9, a feature often cited as one compelling reason for 
using exponential smoothing. It Is true that 9 will have to be 
stored and updated at the next time period and hence one extra 
operation, that of taking the square root, will be involved. But 
with modern ccnaputers, the time for this extra operation is 
negligible. 

The actual point by point results of the simulation are perhaps 
not too enlightening. The program that was used has been preserved 
and is available for further use. Of more interest is a comparison 
of the two procedures with regard to bias and variance. These 
have been estimated by computing the sample average and sample 
standard deviation (S.D.) of each of 9 and 9 for each parameter 
choice. While many more cases were examined, the results for 
parameter choices furnished by NavSup are sutmarized in Table 1. 
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Parameter 

Pairs 


Average 

o 


Average 

Q 


Bias 

a 


f''* ' — 

Bias 

o 


S «D • 

a; 

a 


S.D. 

a 


Me S e £ o 

a 


T 

M.S.E. 

o 


(50,10) 


10.27 


9.60 


0.27 


•=0.40 


2.83 


1.33 


8.08 


1.93 


(100,25) 


25.67 


23.99 


0.67 


-1.01 


7.08 


3.31 


50.57 


11.98 


(100,50) 


51.33 


47.98 


1.33 


-2.02 


14.15 


6.63 


201.99 


48.04 


(400,40) 


41.07 


38.39 


1.07 


=1.61 


11.32 


5.30 


129.29 


30.78 


(500,50) 


51.38 


48.52 


1.38 


=1.48 


14.15 


7.40 


202.13 


56.95 


(600,60) 


61.60 


57.58 


1.60 


°2.42 


16.98 


7.95 


290.88 


69.06 


(700,70) 


64.99 


58.78 


=5.01 


>11.22 


21.40 


10.56 


483.06 


237.40 


(800,80) 


81.51 


76.15 


1.51 


■=3.85 


19.18 


5.85 


370.15 


49.04 


(900,90) 


87.26 


86.21 


=2.74 


>3.79 


21.09 


5.01 


452.30 


39.46 


(1,000,100) 


101.89 


95.18 


1.89 


-4.82 


23.97 


7.32 


578.13 


76.81 



Table 1. Sample Characteristics 



The table clearly brings out the inefficiency of the smoothing 
technique compared with classical methods of estimatlono Both proce» 
dures are biased. The smoothing procedure tends to overestimate while 
the maximum likelihood procedure more conservatively underestimates 
CT, a fact known from the theory of course. Except for the one 
extreme case (the pair (700 j, 70) where the bias for a is numerically 
as large as 11 for some reason that is not consistent with the 
other results, the bias in both cases is of the same magnitude. There 
is a marked difference in the variance estimate of the two procedures, 
however. The variance for a is significantly higher than that for 
o in every case. Indeed, there was never a single normal case 
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studied in which this was not the case. 



Since the bias is roughly the same for each procedure this means 
that the mean squared error (M.S.E.) of each procedure will largely 
be determined by its variances. Recalling that M.S.E. is given by 
the variance plus the square of the bias we use the averages and 
standard deviations of Table 1 to compute estimates of M.S.E. which 
are recorded in the last two columns. Again the difference in the 
two procedures are quite striking. The M.S.E. for a reaches as 
high as 12 times that for o in one case and even in the pathological 
case cited above (the pair (700,70)) where the bias of a is 

A 

unexplainably high, the M.S.E. of a is still twice that of a. 

In most of the cases, the ratio is about 4. 



A similar analysis was undertaken for Poisson demand with results 
that are not quite as striking. These are summarized in Table 2 for 
three typical cases. 



i Parameter 


a 


Average 

a 


Average 

A 

o 


Bias 

a 


Bias 

A 

a 


S.D. 

CJ 


S.D. 

A 

a 


M.S.E. 

a 


M.S.E. 

A 

a 


.01 


0.1 


.01 


.06 


-.09 


«,04 


.04 


.04 


.0097 


-.0032 


0.1 


0.3 


.17 


.23 


-.13 


8 

• 

O 


.13 


.06 


.0338 


.0085 


1.0 


1.0 


.91 


.95 


-.09 


o 

• 

8 


.28 

1 


.12 i 

1 1 , , ,< 


.0865 


.0169 



Table 2. Poisson Demand 

Again it should be observed that the M.S.E. for o is at least 3 

A 

times that for a in each instance examined. 
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4. EFFECTS ON REORDER LEVELS 



Not satisfied to merely summarize the evidence obtained as in 
the tables of the preceding section, it was deemed advisable to test 
the effects of the forecasting errors on their use in fixing reorder 
levels. Again for reporting purposes, the main concentration was spent 
on the normal case. If demand is normal with mean ^ and standard 
deviation o, then the (theoretical) reorder level would be set at 
+ ka where k is chosen to satisfy a given risk p defined by 

p = P(X > p. + ka) 

Thus, the risk, for present purposes, is the probability that demand 
will exceed the reorder level causing a stockout. For given p (or k) 
the parameter k (or p) can be determined from standard normal tables. 

Of course and cx are unknown so that p can never be 
satisfied exactly. If we were to use and o* for the reorder 

level and it should happen that |i^+kO’^>|i+ka where k has 
been chosen to satisfy a given risk requirement p, then the true 
risk say p^ “ P(X > + ko^) would be something smaller. This 

means that simultaneously overestimating \i> and <J results in over- 
stocking, that is, we could stock less to achieve the required risk 
level. There may also be penalty costs in the way of storage costs 
to consider for such a situation. On the other hand, if + kff* < 
p + kcr then, while we think we are stocking in such a way as to 
achieve a given risk p, in fact the true risk p* would be greater. 
The resulting shortage "cost" paid for such a position could be 
disastrous. Since neither position is particularly favorable, and 
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the true parameter values are unknown, it Is clear that the most 
precise estimates of and o are desirable; precise in the sense 
of minimum fluctuation about the true values. 

Trying to determine the true risk Incurred when exponentially 
smoothed estimates |Ji and a of and <7, respectively, are used 
theoretically requires the joint probability distribution of X, p- 
O'. If finding the distribution of o alone seems difficult, the 
task of finding such a joint distribution looms formidable to say the 
least. Fortunately, in a simulation approach the parameters, and 
hence the true theoretical reorder level for a given risk p, is 
known. For any procedure used to set actual reorder levels it is then 
possible to observe the behavior of repeated applications of such a 
procedure. There are many ways this might be done. One approach is 
to observe the number of times a reorder level falls within so many 
units of the theoretical level. Or, for a given percentage P we 
might ask within how many units of the theoretical reorder levels 
will P percent of the actual reorder levels be found? We have 
done a little of both and summarized the findings for various cases 
in Table 3. 

The table is constructed as follows. First a triple of parameters 
P) is chosen where p is the desired risk. From this k is 
determined from tables and a theoretical reorder level computed. For 
example, if M> = 400, a = 40 and p = .01 then from .01 “ P(3^400+40k) 
we determine k ® 2.326 and theoretical reorder level is 493.04. Then, 
for each parameter triple Qj^ is the percentage (rounded) of 200 
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reorder levels found to be within 1 unit of the theoretical reorder 

A? ^ 

level first for the level based on exponential smoothing |x + ko 
and 9 seconds for a level set using X + kc where o is the maximum 
likelihood estimate of o and X is the average demand. Next, P^q 
represents the number of units about the theoretical reorder level 
within which 50% of the computed reorder levels were found. Again, 
the first column under P^q is for smoothed estimates and the second 
for maximum likelihood estimates. The last column is a similar 



computation for 95% of the reorder levels. 



Parameter 

Triples 





^so 


^95 


(100, 25,. 01) 


4 


32 


■ 15 


2 


35 


13 


(100, 25,. 11) 


6 


48 


11 


1 


24 


6 


(100,25,. 50) 


8 


19 


6 


2 


17 


5 


(400, 40,. 01) 


2 


22 


25 


2 


57 


21 


(400,40, . 11) 


5 


38 


18 


2 


38 


9 


(400,40, .50) 


5 


10 


1© 


3 


27 


7 


(600,60, .01) 


1 


15 


37 


4 


85 


31 


(600,60 , . 11) 


5 


26 


26 


3 


58 


14 


(600,60, .50) 


3 


6 1 


15 


4 


40 


12 


(800, 80,. 01) 


2 


15 


35 


5 


103 


46 


(800, 80,, 11) 


2 


16 


28 


3 


69 


30 


(800, 80,. 50) 


2 


23 


21 


2 


56 


23 


(1000, 100,. 01) 


2 


13 


44 


6 


129 


58 


(1000,100,. 11) 


1 


12 


35 


4 


85 


37 


, (1000, 100,. 50) 


2 


18 


j 27 


3 


70 


30 



Table 3. Percentiles for Reorder Levels 
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Many mor« cases than those reported In the table were examined 
but the results are omitted for the sake of brevity. The results in 
these other cases displayed precisely the same pattern however. 

Indeed; the consistency of the various cases is somewhat startling. 
Once again there is not one single choice of a mean, a standard 
deviation or a risk in which the smoothing procedure does not fall 
significantly short of maximum likelihood techniques. This was to 
be expected of course since the variance of the estimators involved 
determine the amount of flucutation from sample to sample. 

One other interesting observation in Table 3 is the monotonlc 
nature of P^q and ^ 95 * ThuS; while it always takes more units of 
distance (roughly twice as many) to Incltide 957 , of the values using 
smoothing versus maxinmm likelihood; the comparison at the 50% level 
Is even more striking. In addition to requiring roughly ten times as 
many units to cover 50% of the observations for smoothing, the 
number of units for the same percentage using maximum likelihood is 
comfortably small. 

Another way to view the results is to try and estimate the true 
risk Incurred when a given procedure Is followed. One way to do this 
is to imagine that, after 2 C 0 observations, the estimates say M> 
and CT* coincide with the parameters and then compute P( 3 ^p>*+k<J*) 
as a measure of the true risk. With a sample size of 200, the 
asymptotic result of supposing that estimates and parameters are the 
same should be fairly sound. Of course the point is that how sound 
such an estimate is depends upon the precision of the estimate once 
again. In any event in practice this is essentially what is done and 
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200 is a figure quite a bit larger than the. typical number of periods 
for which historical data ate araintained for making such estimates 
within NavSup. 

In Table 4 our results are summarized as follows. Again, certain 
triples of parameters It, o and p are selected. This time, for 
each pair (p., a) down a column, five different risks are associated 
across a row. For each such combination, an actual reorder level is 
determined in two ways. First of all, k is determined and then 

^ «s> ^ fit 

reorder levels fixed at p + kcr and X + ko. After the 200th 
observation = P(X>p*+ko’*) is computed and entered In each cell 
(suitably rounded), first for (p and second for 

(P ) “ (X,o) . For lack of a better name, p is called the actual 
risk. In this way, the actual risk incurred using smoothing may be 
compared with maximtan likelihood procedures. 
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.13 
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.52 


(400,40) 
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.13 
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,30 
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. 60 
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.13 


,11 


.30 


.25 


.60 


,52 


(600,60) 


.01 


.01 


,05 


,05 


,13 


,11 


.30 


,25 


.60 


.52 


(800,80) 


.06 


.01 


.18 


,05 


,31 


.12 


.52 


.26 


.76 


.52 


(900,90) 


.12 


.02 


.20 


,07 


.27 


,13 


.37 


.27 


.50 


.51 


(1000,100) 


.06 


.01 


.18 


,05 


.31 


,12 


,52 


,26 


,77 


.52 



Table 4. Actual Risks 



The results are fairly self=explanatory. The results are even fairly 
consistent in the sense that for most parameter choices the asymptotic 
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results are quite the sante. For large means and variances more 
variation begins to show up. We see that for small values of p, both 
precedures are roughly equivalent. That Is to say^ after 200 
predictions;, the actual risk Is fairly close to the true risk. But 
for large or even moderate values of p It appears that smoothing 
consistently produces a larger actual risk than maximum likelihood. 
Both are larger than the true risk but the maximum likelihood 
procedure is much closer to the true value and Is far more consistent. 
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H + ka 


+ ka 


fit 

X + ka 


.01 


51 


38 


44 


,02 


52 


38 


43 


o 
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52 


40 


41 
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53 


39 


40 


.05 


54 


39 


39 


.06 


52 


38 


37 


.07 


52 


40 


37 


.08 


52 


40 


34 


.09 


52 


40 


35 


o 

o 


51 


41 


34 



Table 5. Percentage Overestimates 
Finally » supposing that it is far more serious to overestimate 
the reorder level (resulting in excess stock for a given risk) than 
to underestimate it^ the data were analyzed to see the percentage of 
time the reorder level was cverestimated in 200 trials. For this 
illustration only one distribution was used, namely, normal with 
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mean 50 and standard deviation 10. Again^ for various risks 
(this time no more than 10%) reorder levels are set using + ka 

ca 

and X + ka • This time we examined a third alternativej> a mixture 
of the two procedures s given by |Ji + ka. The results are summarized 
in Table 5 and the entries are the percentage of times that the given 
method overestimates the theoretical reorder level. 

It is by this time perhaps not surprising to note that smoothing 
consistently overestimates the reorder level. Even the mixed procedure, 
whereby smoothing is used to estimate the mean and o is estimated by 
maximum likelihood, is an improvement over a MAD estimate of a. 

Compared with straight maximum likelihood, it is better for small 
values of p and then begins to fall off as p increases. 

5. CONCLUSIONS AND RECOMMENDATIONS 

In drawing conclusions it would be well to reiterate the tentative 
nature of simulation results. Without the exact probability distri- 
butlon for the exponentially smoothed estimators treated in this 
report“”particularly MAD*no definite statement can be made about the 
efficiency of this method over others. But the simulation results 
display a certain consistency that certainly lends support to the 
claim that maximum likelihood methods would be superior for those 
purposes to which smoothing is presently employed. Certainly the 
evidence is sufficient to warrant further examination and comparison 
with alternatives. In spite of the difficulties involved in the 
derivation of the probability distributions when exponential smoothing 
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is employed as a forecasting tool, it is our recommendation to continue 
the theoretical search for those distributions or at least the moments 
that are involved. In particular it would be helpful to establish 
that the variance of o is bounded below by a positive quantity as 
seems to be indicated by the simulations. If this were true then the 
estimate will always fluctuate about the true o regardless of how 
many time periods are considered. 

Another reason for continuing research in this are* is the fact 
that even the simulation studies were only carried out for the normal 
case (and the isolated Poisson cases) with constant s^iean and variance. 
This is a strong assumption and entails supposing that no trend is 
present. It is highly recommended that a similar analysis in the 
presence of trend be carried out. The trend test and indicator 
presently used by NavSup as outlined in C4] again employs exponential 
smoothing as a basic tool for ar.alyzlng the data. If the lack of 
efficiency indicated in this study carries over to that case as well 
then the consequences should be even more serious. Indeed, there is 
reason to believe that with the methods presently employed, there 
can be a 50% chance of rejecting an almost vertical trend, that is 
a stidden jump in the demand pattern. These matters are certainly 
topical candidates for further research. 

In any or all of the above recommendations, it is suggested that 
real data be used to analyze, unfortunately perhaps in retrospect, 
the reorder system had methods such as maximum likelihood been used. 
Large amounts of such data are avaiJ^ble for a wide variety of 
inventory items and the programs for processing such data have been 
preserved. 
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6. APPENDIX 



This appendix Is devoted to supplying the mathemttlcal details 
for the monotonlc nature of ^ as asserted In Section 3. It should 
be recalled from that section that ^ , expressed as a function of 
the Poisson parameter X Is given by 



A 

a 

In particular, when X = n 
u^ given by, 



2 A 



•X X' 



Cx] 



[X]J 



— has a relative minimum with a value 
a 
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n 



2 /n 



=n n 



n 



n! 



It Is conjectured that this sequence of values is monotone increasing. 
To see this observe that 



u 

n 



'^n+l 




and 



log = (n4%) log + 1 « =(n44^) log(l+^ + 1 . 

n+1 

Recalling the Taylor expansion of log (1+x) (see [6 3 example), 
we have, for n > 1, 

1 1 1 1 

log il+3 = n " TIE + R 
® “ 2n 



where 



|r| 
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Then. 






= + R] + 1 

In 2 n 4n2 3 q2 6h3 2 



1 1 2nrfl 

I2n^ 6n^ 2 



R 



But R < 0 for this case so that 



2n+2 2n+2 1 . ^ 1 1 

— “““ 2 “ 1*1 *~3 ■"“4 

^ 4n-^ 8n^ 



Hence, 



log 



u 



n+1 



1 1 . 1.1 

5 " ^ ■*" — S 

12n 6n 4n 8n 



(2n -2n-3> 
24n 



Now for X i 2, the polynomial 2x^ ■=• 2x - 3 is positive so that 



u u 

log ° < 0 anid hence ° < 1 
“n+1 "n+1 



In other words. 



^ ^ ^ J .1 



n “ 2 s35*** 



which is what was to be verified. 
From [ 7 ] (page 64 ) 5 



,-(n+i)^ , 24(n+%) ^ ^ ^ , -(n-^ 
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we have 



c / j.i.\ 2 , “(n+%) 

Since f (n+5) = —7 (n+%) e 
n n* 

1 

/l < f (n+ 4 ) < /2 e 
TT " TT 

Thus the sequence of relative maxima of ^ converges to 
the claim on page 7 that — then converges to this limit. 



supporting 



- 26 - 



BIBLIOGRAPHY 



[l] Zehna, Peter W., Some Remarks on Exponential Smoothing . 

Naval Postgraduate School, Technical Report No. 72, December 1966. 
C 23 Ashar, V.G. and T.D. Wallace, "A Sampling Study of Minimum 

Absolute Deviation Estimators," Operations Research . Vol. 11, 

No. 5, September=October 1963. 

[33 Brown, R.G., Smoothing, Forecasting and Prediction of Discrete 
Time Scries , Prentice “Hall, Inc., 1963. 

[4] PAR 1 - Application D, Operations 5, 6, Levels Computation for 
Consumables and Repalrables, FMSO, Mechanlcsburg, Pennsylvania, 
1964. 

Cs] Held, A., Statistical Theory with Engineering Applications . 

John Wiley & Sons, 1952. 

Cs3 Taylor, Angus E., Advanced Calculus , Ginn and Company, 1955. 

C?3 Feller, William, An Introduction to Probability Theory and Its 
Applications, Vol. 1, John Wiley & Sons, 1950. 



INITIAL DISTRIBUIION LIST 



No. 



Defense Documentation Center (DDC) 

Cameron Station 
Alexandria, Virginia 22314 
Attn: IRS 

Library 

Naval Postgraduate School 
Monterey, California 93940 

Captain Thomas J. Ingram, III, SC, USN 
Director, Research and Development Division 
(Code SUP 063) 

Naval Supply Systems Command 
Department of the Navy 
Washington, D* C. 20360 

Mr. James W. Prichard 
(Code SUP 04E) 

Naval Supply Systems Command 
Department of the Navy 
Washington, D. C* 20360 

Commander Jack A* White, SC, USN 
Fleet Material Support Office 
(Code 97) 

Naval Supply Depot 
Mechanicsburg, Pennsylvania 17055 

Professor Peter W. Zehna 
Department of Operations Analysis 
(Code 55Ze) 

Naval Postgraduate School 
Monterey, California 93940 

Professor David A. Schrady 
Department of Operations Analysis 
(Code 55So) 

Naval Postgraduate School 
Monterey, California 93940 

Professor Jack R. Borsting 
Department of Operations Analysis 
(Code 55) 

Naval Postgraduate School 
Monterey, California 93940 



of Copies 
20 

2 

5 



2 

5 

25 

1 

1 



1 



Commander R. E« DeWlnter, SC„ USN 
Naval Aviation Supply Office 
700 Robb in Avenue 
Philadelphia, Pennsylvania 

Professor K. T« Wallenius 1 

Department of Mathematics 

Clemson University 

Clemson, South Carolina 29631 

Lieutenant W» E» Daeschner, SC, USN 1 

PERS Allb 

Bureau of Naval Personnel 
Arlington Annex 
Washington, D. C. 20360 

Dean C. E. Menneken 2 

Dean of Research Administration 
(code 023) 

Naval Postgraduate School 
Monterey, California 93940 



UNCLASSIFIED 

SooirUv Classificalion 



DOCUMENT CONTROL DATA - R & D 

5«*cufifK clfiisihcation of title, body ol ttbslntct nnd indexintj finnolfttion munt be e/ifered whrn the overuli rrport tn c inuntHcd) 



I ONIGINATINC A C T I V I T Y f CofpOf A au /hor; 

Naval Postgraduate School 
Monterey, California 



2a. REPORT security CLASSIFICATION 

UNCLASSIFIED 



Zb. GROUP 



J REPORT title 



’•forecasting errors USING MAD" 



4 DESCRIPTIVE notes (Type of report end, inclusive detes) 

Research Report 



S AU ThORiSI (First nmme, middle initiel, lest neme) 

Peter W. Zehna 



e repor t date 

April 1969 


7a. total no. of pages 

33 


7b. NO OF REFS 

7 


6A CONTRACT OR GRANT NO 

NAVSUP RDT&E No. TF015-02-101 

6. P ROJ EC T NO 


9e. ORIGINATOR'S REPORT NUMBER(S) 

NPS 55Ze904lA 


c. 

d. 


9fe. OTHER REPORT NO(SI (Any other numbers thmt mey be ms si gned 
this report) 


10 DISTRIBUTION STATEMENT 

Distribution of this document is unlimited. 




M SU P PL EM EN T A R Y NOTES 


12 SPONSORING MILITARY ACTIVITY 




Research and Development Division 
Naval Supply Systems Comnand 



13 ABSTR AC T 



A Study of the effects of using meen absolute deviation (MAD) 
to estimate variability in setting reorder levels for the inventory 
of a stock item. The method presently employed by NavSup in setting 
such reorder levels involves exponentially smoothed estimates of the 
mean and variance of the demand process. Any error involved in setting 
reorder levels results in a change in the underlying risk which in turn 
can be translated into costs. Such errors for the method of estimation 
presently employed are compared with standard maximuni likelihood 
procedures. By simulating several nomal systems, the smoothing 
technique is found to be inferior to classical methods with no reduction 
in computational difficulties. 



FORM 



* 1473 

I NOV 68 I “ # 

S/N 01 01 -807-681 1 



DD 



(PAGE 1) 



unclassified 



Security Classification 



s-3l40a 



UNCLASSIFIED 



So -urjtv Cla‘;Rif K'ation 



KF. Y VWOMOS 


LINK A 


LINK B 


LINK C 


ROLE 


A T 


ROLF 


W T 


ROLL 


Af T 


INVENTORY THEORY 
EXPONENTIAL SMOOTHING 
FORECASTING 

MEAN ABSOLUTE DEVIATION 
REORDER LEVELS 















) /r.,1473 (BACK) UNCLASSIFIED 



1 ' 1 - • 



1 



Security Classification 



A - • 1 4 C - 



